Generic Dynamic Scaling in Kinetic Roughening 
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We study the dynamic scaling hypothesis in invariant surface growth. We show that the existence 
of power-law scaling of the correlation functions (scale invariance) does not determine a unique 
dynamic scaling form of the correlation functions, which leads to the different anomalous forms of 
scaling recently observed in growth models. We derive all the existing forms of anomalous dynamic 
scaling from a new generic scaling ansatz. The different scaling forms are subclasses of this generic 
scaling ansatz associated with bounds on the roughness exponent values. The existence of a new 
class of anomalous dynamic scaling is predicted and compared with simulations. 



The theory of kinetic roughening deals with the fate of 
surfaces growing in nonequilibrium conditions . In a 
typical situation an initially flat surface grows and rough- 
ens continuously as it is driven by some external noise. 
The noise term can be of thermal origin (like for instance 
fluctuations in the flux of particles in a deposition pro- 
cess), or a quenched disorder (like in the motion of driven 
interfaces through porous media). A rough surface may 
be characterized by the fluctuations of the height around 
its mean value. So, a basic quantity to look at is the global 
interface width, W{L,t) = {[h{x,t) — h]^)^/'^, where the 
overbar denotes average over all a; in a system of size L 
and brackets denote average over different realizations. 
Rough surfaces then correspond to situations in which 
the stationary width W{L, t — > oo) grows with the system 
size. Alternatively, one may calculate other quantities re- 
lated to correlations over a distance I as the height-height 
correlation function, G{1, t) — {[h{x + l,t) — h{x, t)]^), or 
the local width, 1^(1,1) = {{[h{x,t) - {h)if)i)^/'^, where 
(• • denotes an average over x in windows of size I. 

In absence of any characteristic length in the prob- 
lem growth processes are expected to show power-law 
behaviour of the correlation functions in space and time 
and the Family- Vicsek dynamic scaling ansataz |^,0j2j 

w{L,t)^r/^f{Lim). (1) 

ought to hold. The scaling function f{u) behaves as 

, . / if u -Cl 

•'^"^ ^ \ const, if u> 1 ' ^ > 

where a is the roughness exponent and characterizes the 
stationary regime, in which the horizontal correlation 
length ^(i) '-^ t^/^ {z is the so called dynamic exponent) 
has reached a value larger than the system size L. The 
ratio (3 — a/ z is called growth exponent and character- 
izes the short time behavior of the surface. As occurs in 
equilibrium critical phenomena, the corresponding crit- 
ical exponents do not depend on microscopic details of 
the system under investigation. This has made possible 
to divide growth processes into universality classes ac- 
cording to the values of these characteristic exponents 
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A most intringuing feature of some growth models is 
that the above standard scaling of the global width dif- 
fers substancially from the scaling behaviour of the local 
interface fluctuations (measured either by the local width 
or the height-height correlation). More precisely, in some 
growth models the local width (and the height-height cor- 
relation) scales as in Eq.(|l|), i.e. w{l,t) — f aQ- / £,{t)) ^ 
but with the anomalous scaling function 

I const if ' 

where the new independent exponent a/oc is called the 
local roughness exponent. This is what has been called 
anomalous roughening in the literature, and has been 
found to occur in many growth models P-p^ as well 
as experiments |l^. Moreover ,it has recently been 
shown 1^,0 that anomalous roughening can take two 
different forms. On the one hand, there are super-rough 
processes, i.e. a > 1, for which always aioc = 1- On the 
other hand, there are intrinsically anomalous roughened 
surfaces, for which a/oc < 1 and a can actually be any 
a > aioc- 

Anomalous scaling implies that one more independent 
exponent, aioc, may be needed in order to asses the uni- 
versality class of the particular system under study. In 
other words, some growth models may have exactly the 
same a and z values seemingly indicating that they be- 
long to the same universality class. However, they may 
have different values of a/oc showing that they actually 
belong to distict classes of growth. As for the experi- 
ments, only the local roughness exponent is measurable 
by direct methods, since the system size remains nor- 
mally fixed. Fracture experiments [p^ in systems of vary- 
ing sizes have succeded in measuring both the local and 
global roughness exponents in good agreement with the 
scaling picture described above. 

In this Letter we introduce a new anomalous dynamics 
in kinetic roughening. We show that, by adopting more 
general forms of the scaling functions involved, a generic 
theory of dynamic scaling can be constructed. Our the- 
ory incorporates all the different forms that dynamic scal- 
ing can take, namely Family- Vicsek, super-rough and 
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intrinsic, as subclasses and predicts the existence of a 
new class of growth models with novel anomalous scaling 
properties. Simulations of the Sneppen model (rule A) 
[ p^ of self-organized depinning (and other related mod- 
els) are presented as examples of the new dynamics. 

Firstly, let us consider the Fourier transform of the 
height of the surface in a system of size L, which is given 
by h{k, t) = SajI^C^' ~ ^(^)] exp(jfcx), where the 

spatial average of the height has been substracted. The 
scaling behaviour of the surface can now be investigated 
by calculating the structure factor or power spectrum 

S{k,t) ^ (h{k,t)h{^k,t)), (4) 

which is related to the height-height correlation function 
G{l,t) defined above by 

G{l,t)^^ J2 [l~cos{kl)]S{k,t) 

2-K I L<k<ix I a 

a / — [l-cos(fcO]5(fc,t), (5) 

J27r/L 27r 

where a is the lattice spacing and L is the system size. 

In order to explore the most general form that kinetic 
roughening can take, we study the scaling behaviour of 
surfaces satisfying what we will call a generic dynamic 
scaling form of the correlation functions. We will con- 
sider that a growing surface satisfies a generic dynamic 
scaling when there exists a correlation lenght ^(t), i.e. 
the distance over which correlations have propagated up 
to time t, and ^{t) ~ t^^^, being z the dynamic ex- 
ponent. If no characteristic scale exists but ^ and the 
system size L, then power-law behaviour in space and 
time is expected and the growth saturates when ^ ^ L 
and the correlations (and from Eq.(|^) also the structure 
factor) become time- independent. The global roughness 
exponent a can now be calculated in this regime from 
G{1 ^ L,t :^ L^) L^a (oj. VF(L,< > L^) - L"). In 
general, as we will see below, the scaling function that 
enters the dynamic scaling of the local width (or the 
height-height correlation) takes different forms depend- 
ing on further restrictions and/or bounds for the rough- 
ness exponent values. These kind of restrictions are very 
often assumed and not valid for every growth model. For 
instance, only if the surface were .self-affine saturation of 
the correlation function G{1, t) would also occur for inter- 
mediate scales I at times t l'^ and with the very same 
roughness exponent. However, the latter does not hold 
when anomalous roughening takes place as can be seen 
from the scaling of the local width in Eq.(^. 

Our aim here is to investigate all the possible forms 
that the scaling functions can exhibit when solely the ex- 
istence of generic scaling is assumed. So, if the roughen- 
ing process under consideration shows generic dynamic 
scaling (in the sense above explained), and no further 



assumptions (like for instance surface self-affinity or im- 
plicit bounds for the exponents values) are imposed, then 
we propose that the structure factor is given by 

where the scaling function has the general form 

and the exponent as is what we will call the spectral 
roughness exponent. This scaling ansatz is a natural 
generalization of the scaling proposed for the structure 
factor in Ref. ||l^,0 for anomalous scaling. 

In the case of the global width, one can make use of 

W'{L,t)^-Y,S{k,t) = J ^Sik,t), (8) 

to prove easily that the global width scales as in Eqs.(|^) 
and (H), independently of the value of the exponents a 
and as- 

However, the scaling of the local width is much more 
involved. The existence of a generic scaling behaviour 
like (|^) for the structure factor always leads to a dynamic 
scaling behaviour, 

w{l,t)^ ./G{i:T)=t^g{l/0 (9) 

of the height-height correlation (and local width) , but the 
corresponding scaling function g{u) is not unique. When 
substituting Eqs.(|^) and (||) into (||), one can see that the 
various limits involved (a —> 0, ^(t)/L oo and L —^ oo) 
do not conmutc [ p^|pT[ . This results in a different scaling 
behaviour of g{u) depending on the value of the exponent 
as. 

Let us now summarize how all scaling behaviours re- 
ported in the literature are obained from the generic dy- 
namic scaling ansatz (Q). We shall also show how a new 
roughening dynamics naturally appears in this scaling 
theory. Two major cases can be distinguised, namely 
as < 1 and > 1. On the one hand, for as < 1 the 
integral in Eq. (||) has already been computed [|l6|,0 and 
one gets 

, . / U°'' if 7i < 1 , . 

^^-<i(")~|const ff u»l- (10) 

So, the corresponding scaling function is 5q^<i ^ /a and 
as — aioc, i-B- the intrinsic anomalous scaling function 
in Eq.(^). Moreover, in this case the interface would sat- 
isfy a Family- Vicsek scaling (for the local as well as the 
global width) only if a = as were satisfied for the par- 
ticular growth model under study. Thus, the standard 
Family- Vicsek scaling turns out to be one of the possible 
scaling forms compatible with generic scaling invariant 
growth, but not the only one. 
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On the other hand, a new anomalous dynamics shows 
up for growth models in which > 1. In this case, 
one finds that, in the thermodynamic limit L oo, the 
integral Eq.(^ has a divergence coming from the lower 
integration limit. To avoid the divergence one has to 
compute the integral keeping L fixed. We then obtain 
the scaling function 



9as>i{u) 



u if 7i <C 1 
const if u ^ 1 



(11) 



So that in this case one always gets aioc = 1 for any 
Us > 1. Thus, for growth models in which a — as one 
recovers the super-rough scaling behaviour p6| , p^ . 

However, it is worth noting that neither the spectral 
exponent as or the global exponent a are fixed by the 
scaling in Eqs(|^) and ( ITTl) and, in principle, they could 
be different. Therefore, growth models in which ag > 1 
but a ^ as could also be possible and represent a new 
type of dynamics with anomalous scaling. The main fea- 
ture of this new type of anomalous roughening is that it 
can be detected only by determining the scaling of the 
structure factor. Whenever such a scaling takes place 
in the problem under investigation the new exponent as 
will only show up when analyzing the scaling behaviour 
of S(k,t) and will not be detectable in either W{L,t), 
w{l,t) or G{l,t). In fact, as we have shown, the station- 
ary regime of a surface exhibiting this kind of anoma- 
lous scaling will be characterized by W{L) ^ and 
1^(1, L) ^ G{1, L) ^ IL""^^, however, the structure 
factor scales as S{k,L) ~ k-i^a^+i) i2{a-a,) ^j^ere the 
spectral roughness exponent as is a new and indepen- 
dent exponent. We can summarize our analitycal results 
as follows 



if as < 1 ^ aioc = as 
if as > I ^ aioc = 1 



as = a ^ Family — Vicsek 
cks ^ a => Intrinsic 

= a => Super — rough 
as ^ a ^ New class 
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In the following we present simulations of a one- 
dimensional growth model that is a nice example of the 
new dynamics. We have performed numerical simula- 
tions of the Sneppen model of self-organized depinning 
(model A) We have found that this model exhibits 
anomalous roughening the type described by Eq.(|^) for 
as > 1 and as ^ a. In this model the height of the 
interface h{i,t) is taken to be an integer defined on a 
one-dimensional discrete substrate i — 1, ■ ■ ■ , L. A ran- 
dom pinning force 77(1, h) is associated with each lattice 
site [i,h{i)]. The quenched disorder r]{i,h) is uniformly 
distributed in [0, 1] and uncorrelated. The growth algo- 
rithm is then as follows. At every time step t, the site ig 
with the smallest pinning force is chosen and its height 
h{io, t) is updated /i(io, t + l) = h{io, t) + l provided that 



the conditions \h{io,t) — h{io±l,t)\ < 2 are satisfied. Pe- 
riodic boundary conditions are assumed. We have stud- 
ied the behaviour of the model in systems of different 
sizes from L = 2^ up to L = 2^^. From calculations of 
the saturated global width W{L) for various system sizes 
we find a global roughness exponent a — 1.000 ± 0.005 
in agreement with previous simulations p^ . We have 
checked that the scaling of the global width is given by 
Eq.(|l]) with a scaling function like (||). Also in agree- 
ment with previous work ||l8| we find that the time ex- 
ponent a/z = 0.95 ± 0.05. The local width w{l,t) scales 
as w{l, t) = t'^^'^gil/C) where the scaling function is given 
by Eq.([ll|), and also aioc = 1- 

From these simulation results one could conclude that 
the behaviour of the Sneppen growth model is rather 
trivial and that the exponents a = aioc = z — 1 describe 
its scaling properties. Quite the opposite, this model ex- 
hibits no trivial features that can be noticed when the 
structure factor is calculated. In Figure 1 we show our 
numerical results for the structure factor S{k, t) in a sys- 
tem of size L = 2048. Note that in Figure 1 the curves 
S{k,t) for different times are shifted downwards reflect- 
ing that a < as- This contrasts with the case of intrin- 
sic anomalous roughening ||l^,|l^ where as = aioc and 
c^ioc < 1- The slope of the continuous line is —3.7 and 
indicates that a new exponent as = 1.35 enters the scal- 
ing. This can be better appreciated by the data collapse 
shown in Figure 2, where one can observe that, instead 
of being constant, the scaling function s{u) has a nega- 
tive exponent u^"-^ for u ^ 1. The exponents used for 
the data collapse are a = 1, 2; = 1 and the scaling func- 
tion obtained is in excellent agreement with Eq.(^ and 
a spectral exponent as = 1.35 ± 0.03. 

The interface in the Sneppen model A is formed by 
facets with constant slope ±1 jis). The value of the 
exponents a — aioc — 1 and a^ = 1.35 is related to 
the faceted form of the interface at saturation. It is 
easy to understand how the anomalous spectral rough- 
ness exponent appears due to the faceted form of the 
interface. For the simpler (and trivial) case of a faceted 
interface formed by a finite number of identical segments, 
N, of constant slope, ztm, one can show analitically that 
the global width W{L) m^L^/N^, and the height- 
height correlation function G{1) ^ Pni^ — Nni^l^/L, 
which leads to a = aioc — 1, while the spectrum 
S{k, L) ^ k~'^L~^ as fc — > 0. A simple comparison with 
the anomalous scaling form for the stationary spectrum 
S{k, L) - /c-(2a=+i)7,2(a--a,) ^^g^^g ^ ^ ^ Actually, 

the facets occuring in the Sneppen model are not formed 
by identical segments, but rather follow a random distri- 
bution j20j, which leads to a spectral exponent different 
from the trivial case. 

In summary, we have presented a generic theory of 
scaling for invariant surface growth. We have shown 
that the existence of power-law scaling of the correla- 
tion functions (scale invariance) does not determine a 
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unique form of the scaling functions involved. This leads 
to the different dynamic scaling forms recently observed 
in growth models P-pX|] and experiments [Tl|-[l5|] exhibit- 
ing anomalous roughening. In particular interface scale 
invariance does not necessarily imply Family- Vicsek dy- 
namic scaling. We have derived all the types of scal- 
ing (Family- Vicsek, super-rough and intrinsic anoma- 
lous) from a unique scaling ansatz, which is formulated 
in the Fourier space. The different types of scaling are 
subclasses of our generic scaling ansatz associated with 
bounds on the values that the new spectral roughness 
exponent may take. This generalization has allowed 
us to predict the existence of a new kind of anomalous 
scaling with interesting features. Simulations of a model 
for self-organized interface depinning have been shown 
to be in excellent agreement with the new anomalous 
dynamics. It has recently been shown [l9|] that anoma- 
lous roughening stems from a non trivial dynamics of the 
mean local slopes {{Vh)^). In contrast, the new anoma- 
lous dynamics can be pinned down to growth models in 
which the stationary state consits of faceted interfaces. 
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FIG. 1. Structure factor of the Sneppen model for interface 
depinning at different times. The continuous straight line is 
a guide to the eye and has a slope -3.7. Note the anomalous 
downwards shift of the curves for increasing times. 
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FIG. 2. Data collapse of the graphs in Fig. 1. The expo- 
nents used for the collapse are a = 1.0 and z = 1.0. The 
straight lines have slopes —0.7 (solid) and 3.0 (dashed) and 
are a guide to the eye. The scaling function is given by Eq. 
(7) with a spectral roughness exponent Us = 1.35 ±0.03. The 
deviations from the scaling for large values of the argument 
kt^''^ are due to the finite lattice spacing. 
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